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We present a proposal for entanglement purification of the continuous-variable 
quantum state of two propagating optical fields. The scheme is based on each 
field interacting with a local node -atomic ensemble- whose internal collective 
excitation plays the role of an ancillary continuous variable resource. Entanglement 
purification is achieved by a dichotomic measurement, representing the required 
non-Gaussian element, which consists in detecting the presence or absence of 
collective excitations in the atomic ensemble. This scheme can be extended to 
networks, where the nodes are single trapped atoms, and constitutes an important 
building block for the implementation of a continuous-variable quantum repeater. 
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1. Introduction 

The property of entanglement lies at the heart of quantum information processing and quan- 
tum mechanics in general [JJ. Continuous variable (CV) entangled states are very attractive 
in this regard due to the relative easiness in their generation and manipulation, as compared 
to their discrete counterpart [2]. However, the degree of entanglement between two distant 
sites of a quantum network usually decreases exponentially with the length of a connecting 
channel, calling for the implementation of some entanglement purification procedures [3]. 
This involves purifying mixed entangled states, such as those that are obtained when the two 
halves of the paradigmatic CV entangled state, the two- mode squeezed state [2], have been 
distributed through noisy channels. For this purpose, we recall that a procedure involving 
solely Gaussian local operations and classical communications does not lead to increase of 
the initial entanglement [3], and that non-Gaussian local operations must necessarily take 
part of any entanglement purification protocol [5]. 

In this article we propose a scheme for increasing the CV entanglement of two traveling 
optical field: each beam is sent to a local node where it interacts in an appropriate way 
with atoms. In the limit in which the photonic field is coherent over a large number of 
atoms, a collective ensemble excitation is created. This excitation can be described by 
a continuous variable, a bosonic mode, and its interaction with the input optical mode 
realizes the first part of a purification protocol. The second important part is the non- 
Gaussian measurement, implemented through a dichotomic measure of presence or absence 
of atomic excitations in the ensemble. An alternative realization for the local node can 
be a single atom tightly confined by a harmonic potential, where the bosonic mode are 
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the center-of-mass excitations, which are interfaced with light by means of the mechanical 
effects of atom-photon interactions [6l-HQ| . In this case, the dichotomic measure consists 
in detecting whether or not the atom is in the ground state of the oscillator [TT]. We 
illustrate the efficiency of the protocol for the case when the optical modes are initially 
prepared in a two-mode thermal squeezed state, and show that the protocol is able to 
increase entanglement of the quantum state of the fields at the output with respect to the 
entanglement they possessed at the input. 

We remark that, apart for quantum communication, the possibility to increase the avail- 
able entanglement is helpful in various scenarios, such as for instance in quantum metrology. 
Here, entangled states allow to achieve the ultimate sensitivity allowed by the Heisenberg 
limit in phase-sensitive measurements, such as those employed in atomic clocks, and in weak 
force detection [12]. 

The paper is organized as follows. In Section II, we introduce the theoretical model and 
discuss the basic assumptions. In Section III we use the characteristic function for the state 
of the system in order to evaluate the fidelity of the scheme. In Section IV we discuss the 
results, i.e., the amount of entanglement of the two optical fields at the output of the whole 
process, measured in terms of the teleportation fidelity, for all the possible conditional 
states. Concluding remarks and discussions are presented in Section V. The appendices 
report mathematical details, which complement the theory presented in Sec. Ill and IV. 
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2. The model 

The scheme for achieving CV entanglement purification of two optical modes is described in 
Fig. [TJ Two entangled light fields, with photon annihilation operators a\ and 02, propagates 
towards distant locations, where each field impinges onto an atomic ensemble. Each ensemble 
is homogeneously broadened and is constituted by a large number of identical atoms, whose 
relevant electronic transitions are composed by two stable states \g) and \e) coupled to a 
common excited state and forming a A configuration. The impinging field is far-off- 
resonantly coupled to the transition \g) — > \i) (|e) — > \i)) with strength g, while a classical 
laser drives the transition |e) — > \i) (\g) — > \i)) with coupling strength £1 (the two possible 
excitation models implement different dynamics, as discussed below). The parameters are 
chosen so that the excited state is effectively empty. The resulting effective interaction is 
a coherent two-photon coupling between the stable states \g) and |e), while spontaneous 
emission processes are negligible. In the following treatment we assume that density and 
size of the ensembles are essentially equal, so that the laser and input-field parameters can 
be taken to be the same for both atomic ensembles. 

2. A. Atom-Photon Interaction 

In this paper the interaction between atoms and light is assumed to be Hamiltonian. Denot- 
ing by p the density matrix of optical fields and atoms, its evolution during the interaction 
is described by the von-Neumann equation 



where the specific form of H may take two forms (parametric or beam-splitter interaction) 
and is given below. For the schemes discussed in what follows, it is assumed that the atoms 
are all initially prepared in the internal state \g). 

2.A.I. The Parametric Interaction 

We first consider the configuration shown in Fig. QJa). Here, the input field couples to 
the atomic transition |e) —> \i), while the laser drives the transition \g) — > \i). In the 
regime in which the excited state is far-off resonance, it can be adiabatically eliminated 
and the dynamics of quantum field and atomic excitations can be described by the effective 
Hamiltonian [T91I20] 

^=E^V^gaJ + H.c.) (2) 

1=1,2 ^ ' 

where the index 1 = 1,2 denotes the input optical mode and the atomic ensemble it impinges 
on, and Seg the collective spin operator for the corresponding ensemble. Here, Ni is the atom 
number in ensemble I and we denoted by A the detuning between the optical fields and the 
excited state. In Eq.((2J), states \g) and |e) are resonantly coupled (ac-Stark shifts have been 
neglected as they can be compensated by suitably tuning the laser frequency). 

In the Holstein-Primakoff representation |22| . the collective atomic operators (identified 
as angular momentum operators) can be associated with bosonic creation and annihilation 
operators via the relation Sgl = {N\ — blfy) 1 / 2 bi/\fN~i, where b\ denotes the bosonic anni- 
hilation operator satisfying the canonical commutation relations [&Zj£>t] = &ik- I n the low 
excitation limit, where the number of atoms transferred to state |e) is small compared to the 
total number of atoms, the collective atomic operators can be approximated by Sgl ~ bi, 
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corresponding to taking the lowest order term in the expansion in powers of {b\bi)/Ni. By 
suitably choosing the phase of the coupling strengths g and Oj, the Hamiltonian ([2]) can be 
reduced to the form 

H = ^2 ih^i{a\b\ - aih), (3) 

1=1,2 

where [ii = y/Ni\gQi\/ A. Therefore, when the number of excited atoms is much smaller 
than the total number of atoms, the collective internal atomic excitation can be treated as 
continuous variable and the light-atom interaction for the setup of Fig. DJa) is that of a 
parametric amplifier, in which a photon and a quantum of the atomic collective excitation 
are simultaneously created (annihilated). We shall assume from now on that Hi = H2 = H 
(which could be realized by suitably adjusting the intensity of the two driving fields). 

2. A. 2. The Beam Splitter Interaction 

We now consider the configuration shown in Fig. UXb). The input optical field drives the 
\g) — > \i) transition of the ensemble, while the laser couples to the transition |e) — > \i). Adia- 
batic elimination of the excited state gives a beam splitter form of the effective Hamiltonian, 
which reads 

H= ^ mu^ajk - ai b\), (4) 
1=1,2 

where v\ = y/Ni\gQi\/ A and we made a suitable choice of the phase of the coupling strengths. 
Here the generation of a photon is accompanied by the absorption of an atomic collective 
excitation. We assume also in this case that 1^1 = ^2 = v- 
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2. A. 3. Discussion 

The two effective Hamiltonians introduced in this section determine the interaction be- 
tween fields and atoms discussed in this paper. The assumption of Hamiltonian dynamics 
is of course correct when incoherent processes can be neglected during the protocol. On the 
one hand, as we already remarked, only stable atomic states are involved in the dynamics, 
while radiative decay can be neglected due to off-resonant excitation. Moreover, it is also 
reasonable to neglect dephasing of the two metastable states which may arise from collisions 
between the atoms of the gas, since for sufficiently dilute media the relevant coherence be- 
tween the levels \g) and |e) can be maintained for a sufficiently long time [131417] . One could 
avoid the effects of atomic collisions by choosing cold atoms, very dilute gases, or defects in 
solid matrices as rare-earth-doped solids, the latter suffering however from inhomogeneous 
broadening |18j . In such a case one has a smaller optical depth which, in our model, means 
smaller coupling constants /U and v, imposing longer interaction times. A further possibil- 
ity is to employ nodes based on single atoms in traps, using the center-of-mass motion as 
bosonic mode. In this case large optical depth can be achieved by means of optical ele- 
ments, such as a fiber or a cavity. In this context we note that the dynamics described by 
Hamiltonian (|3j) and @ can be derived for this specific setup, where a is here the cavity 
mode operator and b the operator annihilating an excitation of the center-of-mass motion, 
see Refs. [SHE]. 

For the case of atomic ensembles, an implicit but important physical assumption of our 
description is perfect spatial mode matching between each optical mode impinging on the 
atomic ensemble and the quantum optical modes naturally excited by the driven atomic 
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medium in the A configuration. Assuming perfect mode matching means that the output 
field mode emerging from the atomic ensemble after the interaction (and the measurement) 
is exactly the same optical mode entering the ensemble and therefore they are described by 
the same bosonic operators dj. In general, however, the spatial properties of the emerging 
field are determined by the geometric configuration of the gas cell and of the driving clas- 
sical beams |191l20j and generally do not perfectly match with those of the input beams. 
Such limitation could be overcome with suitable optics techniques. Significant results with 
entangled single photon beams have been recently achieved [21j . 

2.B. The measurement 

After the interaction with the atomic ensemble, the second step leading to entanglement 
purification is a non-Gaussian operation. The easiest way to implement it is to measure the 
state of the ancillary modes (the atomic collective modes), checking if a bosonic mode is in 
its vacuum state or not. This consists in measuring whether each ensemble possesses or not 
spin excitations. In the case of an optical mode this projective measurement is performed 
by single-photon on-off detectors which check if there is at least one photon in the mode. In 
an atomic system this measurement can be implemented by suitably adapting the electron 
shelving or "quantum jump" technique [25] to this atomic ensemble situation. Here, one 
applies a laser pulse, coupling state |e) with a fourth atomic level. If no fluorescence photon 
is observed it implies that state |e) is not populated, i.e., all the population resides in 
\g) (outcome 0). If scattered photons are observed, some atomic excitations are present 
(outcome 1). The detection procedure hence requires to distinguish no excitations from the 
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presence of at least one excitation. For the atomic ensemble case, this requires measurement 
procedures like the one discussed in Refs. j26j[27]. When the node is a trapped ion, the 
procedure can be implemented by means of transferring the population of the ground state of 
the oscillator from the initial electronic state into a hide qubit and measuring the resonance 
fluorescence of the electronic state in a closed optical transition [11]. The possible outcomes, 
(1) correspond to the projectors |0) b (0| (lb — |0) b (0|). The projectors, corresponding to 
the measurement result, are described by the POVM elements 

{^=0,^=1} = {|0) 6l (0|, I 6l -|0) 6i (0|}, (5a) 
{E X2=0 , E X2=1 } = {|0} 62 (0|, I b2 -|0} 62 (0|}, (5b) 

where \j is the identity operator for the collective spin in ensemble j and |0) 6 (0| the 
projector into the ground state (absence of excitations), which here corresponds to all atoms 
in state \g). The implementation of this step represents the non-Gaussian operation which 
may allow for purifying entanglement. Denoting by x = (xi, X2) the generic outcome, the 
resulting density operator for modes a\ and 02 conditioned to the outcome x can be written 
as 

Px = —Ttbite&Ex}, (6) 

Px 

where p x = Tr aia2 {p x } is the probability of the measurement outcome and Tr?,^ represents 
the trace over the atomic degrees of freedom. The (conditional) density matrix p x describes 
now the state of the fields. The unconditional state, obtained without selecting the outcome, 
is p' = YlxPxPx- In this paper we will compare the entanglement of the fields described 
by the conditional state p x with the initial state of the fields, checking whether the initial 
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entanglement has been increased by the procedure. 

3. Time Evolution of the System and non-Gaussian Measurement 

We now calculate the time evolution of the coupled atomic-optical modes and then determine 
the various conditional states after the measurement. We choose to work with the normally- 
ordered characteristic function which is defined as I23U24I 

(7) 



X (ff, t) = Tr L(t) n e*"J e-"? ff < j , 



with a = (01,02,61,62) representing the set of bosonic operators, and fj = (oji, a%, j3i, fa) 
with ctj,/3j complex numbers. 

We shall consider a specific initial state of the optical fields, namely, a two-mode squeezed 
thermal state, while all the atoms are in the state \g), or equivalently in the ground state 
of the ensemble collective excitations. This yields the following form of the corresponding 
characteristic function 

X(rf,0) =exp|- { - ^ ^ 2 + n tb ^j (\ai\ 2 + |a 2 | 2 )| exp | - (aia 2 + Q*«2)| ; ( 8 ) 

with nth number of thermal photons and A the two-mode squeezing parameter. For = 
the state for the input fields becomes a standard two-mode squeezed state. The equation 
of motion for the characteristic function is obtained by taking the Fourier transform of the 
von Neumann equation ([T|), see [24], and for the parametric and beam-splitter case it takes 
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the form 



dXpar(ff,t) 

at 



H ai/3i + al(3l + a 2 /32 + - «i 



0# 








dxBs(f),t) 
dt 




Xpar(ff,t). 



XBs(f),t). 



(10) 



(9) 



where we used the Hamiltonian in Eq. (|3|) and Eq. @, respectively. The solutions of these 
equations are presented in Appendices [A] and [B] and provide a full description of the state of 
the whole system of four modes, an optical and an atomic excitation mode at each separated 
node, after an interaction time t. Notice that the duration of the light-ensemble interaction 
is controlled by turning on and off the two driving laser fields. 

The state of optical fields and atomic excitations, which evolves from the two-mode 
squeezed thermal state, Eq. (|8j) according to Eq. © or Eq. (fTU|) . is a Gaussian state of 
the four modes, since the initial state is Gaussian and the coherent evolution preserves 
the Gaussian nature. In order to increase the entanglement of the optical modes, hence, 
a non-Gaussian element must be introduced in the manipulation of the relevant optical 
modes [HE]. For the measurement described by the POVM in Eq. ©, where the outcome 
is either positive (presence of atomic excitations in the ensembles after the interaction) or 
negative (all atoms in the initial state), only the positive outcomes lead to non-Gaussian 
states of the optical modes. We will then focus on these results in order to check whether 
entanglement has increased with respect to the input state. 

We proceed by determining the characteristic function of the optical modes, resulting 



11 



from the measurement. In order to to so, we first need the density matrix of the full system 
in terms of the characteristic function, 



d 2 ai d 2 a 2 d 2 (3i d 2 /3 : 



7T 7T 7T 7T 



2 

X{ai,a 2 ,f3i,f3 2 ) x 



x J D t (a 1 ) J D t (a 2 ) J D t (/3i) J D t (/3 2 )e ( - |Ql|2 - |Q2|2 - |/3l|2 - |/32|2)/2 , (11) 

where D(ctj) = exp ja^a] — a*aj j denotes the displacement operator [21]. We use this 
density matrix in order to determine the four characteristic functions, corresponding to the 
four conditional states in Eq. ([6j). The latter can be written as 

Xoo(ai,a2) = — I(ai,a 2 ; 1, 1), (12a) 
Poo 

Xoi(ai,"2) = — [l(ai,a 2 ;l,0)-/(ai,a 2 ;l,l)], (12b) 
Poi 

Xio(a;i,a 2 ) = — [/(ai,a 2 ;0, 1) - 7(ai,a 2 ; 1, 1)] , (12c) 
Pio 

Xn(ai,a 2 ) = — [l(ai,a 2 ;0,0) - 7(ai,a 2 ; 1,0) 
Pn 

-/(oi,a 2 ;0,l)+7(oi,a 2j l,l)] (12d) 
where the probabilities p xlX2 can be obtained from the normalization condition 

Xxix 2 (ai = 0,a 2 = 0) = 1 . (13) 

In Eqs. (|12p we have used 

I(a 1 ,a 2 ;u,v)= [ d(u^i)d(^ 2 )x(ai, a 2 , ft, (3 2 ) e -^ 2 -^ 2 , (14) 
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where the measures are defined as follows 



d 2 [3 



1 



d(u/3x) 



< 



(15a) 



u 







V 



1 



< 



(15b) 



ir5 2 {fo)d 2 fo 



v 







The expressions for the integrals appearing in Eqs. ([12)1 and the accompanying probabilities 
are given in Appendices O and [Dj 

3. A. Evaluating the Entanglement: Fidelity of Teleportation 

We remark that Eq. (|14h is a Gaussian, and correspondingly the characteristic function 
for the outcome (0,0), Eq. (112a ). is a Gaussian state. The characteristic functions for the 
other outcomes, in which at least one excitation is detected in one ensemble, are clearly 
non-Gaussian and could hence possibly be states where the initial entanglement has been 
purified. We note that these conditioned, non-Gaussian states are not pure. For this kind 
of states the available entanglement measures are rather cumbersome. In order to quan- 
tify the entanglement of these states we adopt the teleportation fidelity [28] as operational 
measure. The teleportation fidelity quantifies entanglement in terms of resource for tele- 
portation. Here, the conditional output state is assumed to be initially shared by Alice 
and Bob for the standard CV-teleportation protocol [29H3T] and the average fidelity F for 
teleporting a coherent state is calculated. In terms of characteristic functions, when these 
are symmetrically ordered, the teleportation fidelity takes the form |31j 




(16) 
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where $ ch is the characteristic function for the shared bipartite channel, the optical 
fields at the output in our case, and <E> m (£) = exp(— 1£| 2 /2) the characteristic function 
of the coherent state to be teleported □. Using the relation between symmetrically- and 
normal-ordered characteristic functions |32|. and 



(17) 



for the shared bipartite channel, we arrive at the expression 



(18) 



Inserting the explicit expressions of the normally ordered characteristic function of the four 
conditional states in the right-hand-side of Eq. Q18|) , after performing the integral we obtain 
the four conditional teleportation fidelities, 



F 00 = — x 



Poo (fl + l) 2 -fl? 2 2 + 2A-2A 12 -2^ 



(C+£>) 2 



-Bvz 



Ft 



01 



F 
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Poi 



B + l 2 + 2A-2A 12 -^f 

1 1 



{B + l?~Bl 2 2 + 2A-2A 12 -2^±^- 2 



'li 



1 



2 + 2A-2A 12 B + l 2 + 2A-2A 12 -£+g)l 



+- 



B + iy-B 2 12 2 + 2A-2A 



o (C+D? 
12 A B+l-B 12 



(19a) 



(19b) 



(19c) 



where the coefficients are given in Appendix A and B for the case of parametric and 
beam-splitter interaction, respectively (see also Appendices C and D for their final explicit 
expressions). CV-entanglement purification is achieved if one of these conditional fidelity 



Since in the expression of the fidelity the characteristic function of the coherent state appears with the 
modulus squared, a possible displacement of a coherent state does not matter. 
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is larger than the teleportation fidelity, evaluated when the bipartite channel is the initial 
state ([E|), and given by 

2l-A + n m (l-A 2 )' 1 Uj 

We will denote it by initial teleportation fidelity. 

4. Results 

We evaluate the teleportation fidelity for the parametric and the beam-splitter interaction 
between light and atoms, and study it as a function of the initial two-mode squeezing, 
characterized by the squeezing parameter A, and by the number of thermal photon in the 
initial state nth- We first consider the case of the parametric interaction. Figure [2] displays the 
four conditional fidelities versus the interaction time t, where the initial teleportation fidelity 
is shown for comparison. One clearly observes that none of the four conditional fidelities 
is larger than the initial one, implying that this protocol does not achieve entanglement 
purification when the interaction between light and atoms is parametric-like. This fact 
can be intuitively explained in terms of monogamy of entanglement [33]: at each node 
the parametric interaction generates entanglement between the incident optical mode and 
the atomic collective mode, and this latter entanglement is detrimental for the all-optical 
entanglement we aim at increasing. For the beam-splitter case, an analysis of the equations 
for the teleportation fidelities, Eqs. (|19p . shows that Fqq, Fq\ and -Fio are always below i^ n it- 
On the contrary, Fn, the fidelity corresponding to the measurement outcome 11, in which 
excitations are detected in both ensembles after the interaction can be greater than Fi nit 
depending on both t and A (besides nth)- In particular there exist an optimal interaction 
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time at which F±\ is maximal. It depends on A and nth, an d it is always close to ixjv. 
Figure [3] displays the teleportation fidelities for the case of beam-splitter interaction as a 
function of A and nth, f° r optimal time. One clearly observes that entanglement purification 
is achieved for Fn when the initial value of A ~ 0.5 and the number of thermal photons 
is small. As n t h is increased, the range of values of A where entanglement purification is 
obtained decreases, till it vanishes. These results suggest to consider F\\ — Fi n u as figure 
of merit of the entanglement purification process. A better estimate of the efficiency of the 
protocol should take into account the probability of the desired outcome, and can be made 
by means of the protocol efficiency £, which we define as 

pn {F n - F mit ) if (F n - Fmt) > 
£ = < (21) 

if (F n - Finn) < 

The deviation of Fn from Fi n u is here weighted by the corresponding probability pn that 
the event occurs. Entanglement purification is achieved when £ > 0. Figure (U^a) displays 
the protocol efficiency as a function of time and initial two-mode squeezing parameter A for 
n th = 0, showing that £ is greater than zero in a wide region of parameters. In this case, 
for nth = the initial state is pure, and in the corresponding region of parameters where 
£ > the protocol achieves entanglement concentration. Figure H^b) displays the efficiency 
when nth > 0. Here, one sees that the protocol may achieve entanglement purification for 
appreciable intervals of interaction times vt and initial two-mode squeezing A. Clearly, the 
range of nth over which £ remains greater than zero is finite. 

If the roles of atomic and optical modes are exchanged in the measurement part of the 
protocol, i.e., one performs a measurement of presence or absence of photons in the opti- 
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cal modes after the interaction, then the protocol achieves CV-entanglement between the 
atomic ensembles starting from entangled light modes and it hence performs a kind of en- 
tanglement swapping. Moreover, we find that in this way the entanglement of the ensembles 
conditional state, corresponding to the outcome (1,1), is larger than the initial entangle- 
ment of the optical modes even in the case of parametric interaction, and in particular, it 
can be larger than the entanglement that could be generated using an effective two-mode 
squeezing interactions between the ensembles. 

The calculations involved in this case mimic closely the calculations already performed 
in Sec [3l with the only difference of exchanging the roles of modes dj and b{. This results 
in fidelities equivalent to those of Eqs. (|19p . but with A B and A\2 ^ B\2- Figure [5] 
displays the corresponding conditional teleportation fidelities for the parametric-amplifier 
interaction, showing an increase of F\\ with respect to i^nit f° r small values of A and short 
interaction times t, while it vanishes as time and A are increased. 

Figure [6] displays F\\ versus time for small, fixed values of A. The conditional fidelity is 
compared with the fidelity, when the entangled pair is generated by a parametric-amplifier 
Hamiltonian (F param = (1 + fit)/ 2). The increase is small, and vanishes as nth augments. 

5. Concluding Remarks 

A protocol for CV entanglement purification of two optical fields has been proposed. En- 
tanglement purification is achieved by letting first the fields interact with an atomic node 
each, constituted for instance by the collective excitation of an atomic ensemble, and then 
by performing a dichotomic measurement of the presence of absence of atomic excitations 
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on the atoms. We have considered two types of interactions between fields and atoms: 
parametric-like (where field and atomic excitations are simultaneously created or annihi- 
lated) and beam-splitter-like (where the total number of photons and atomic excitations 
is conserved). It turns out that by using the effective parametric amplifier interaction, it 
is impossible to improve the teleportation fidelity. This is an indication of entanglement 
monogamy in CV setting [33J . On the contrary, by using the effective beam-splitter Hamil- 
tonian we have shown that it is possible to enhance the input entanglement, even if the 
initial weakly entangled state is not pure. From the perspective of this effective interaction, 
our protocol shares analogies with the protocol for generation of entangled states through 
photon subtraction presented in Ref. [34, 35j . In particular, we have extended the analysis 
to input thermal states and we have analyzed the performance in terms of a new global 
efficiency parameter which takes into account also the success probability of the protocol. 
As a byproduct of our analysis we have shown that, if the role of the optical and atomic 
excitation mode are exchanged such that one performs the non-Gaussian measurement on 
the optical modes, one swaps entanglement between optical and atomic modes, thereby in- 
creasing the final entanglement of the atomic excitations over the initial one of the optical 
fields. In this case, we find a small improvement of the fidelity even when the interaction is 
of the parametric-amplifier type. 

The present protocol could be extended to the situation, in which the distant nodes are 
single atoms tightly confined by harmonic traps and inside high-Q resonators. In this case, 
the harmonic motion of the atomic center-of-mass plays the role of the ensemble excitation 
in the protocol, and a beam-splitter or parametric interaction with the input, optical field 
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can be tailored by using an external laser field [6Hl0j . Measurement of whether the atom 
is in the trap ground state or not can be made by means of electron shelving techniques as 
used in laser cooling [TT]. An extensive discussion of the experimental requirements for a 
setup of this kind can be found in [9lll0|. 

In general, the use of atomic nodes in place of optical elements for entanglement purifica- 
tion of light fields provides one additional control tools, and our work is a first step in this 
direction. 
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A. Characteristic Function - Parametric Interaction 

Since the initial state is Gaussian and the involved Hamiltonian is bilinear, the generic state 
at time t is again Gaussian and the general solution of Eq. Q can be written as 

X P ar%t) = exp[-A(|a 1 | 2 + |a 2 | 2 )- J B(|/3i| 2 + |/3 2 | 2 ) 
+A 12 (aia 2 + a\a* 2 ) + B l2 + ft ft) 

+C (ai/3a + atifa + 012P* + a^i) 

+D{a l p l + aiPt + a 2 p2 + a* 2 Pl)]. (22) 

The individual coefficients in the exponent solve the set of six differential equations 

A = 2/j.D, (23a) 

B = 2fiD, (23b) 

D = n(l + A + B), (23c) 

A 12 = -2(iC, (23d) 

flia = -2fiC, (23e) 

C = -M-Bi2 + ^l2), (23f) 
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with the only nonvanishing initial conditions A(0) = + n th and ^4i 2 (0) = j^i- The 

solutions read 



[1 + n th {l - X 2 )] cosh 2 fit 



A(t) = i \_ X 2 -- 1 ( 24a ) 

[1 + n th (l - A 2 )l sinh 2 

B(t) = ± ^- ^ (24b) 

1 — A z 

D(() = 1 [i + ^d-^i^^) (24c) 

Oft) = _|E^ (24d ) 

*.(0 = (24e) 

M - (24f) 



B. Characteristic Function - Beam Splitter Interaction 

Repeating the same argument of the preceding Section, the general solution of Eq. ()10p for 
the beam splitter case can be written as 

XBs(v,t) = exp [-A(|ai| 2 + |a 2 | 2 ) - 5(|/3i| 2 + |/3 2 | 2 ) + A 12 { ai a 2 + a\a* 2 ) + B 12 (fcfo + /3*/3*) 
+D (axpt + a\l3 x + a 2 /3 2 * + a* 2 p 2 ) + C («i/3 2 + a\fc + a 2 p x + a£#)] . (25) 

The individual coefficients solve the two sets of differential equations 

A = -2i>D, (26a) 
B = 2uD, (26b) 
t) = v(A-B), (26c) 
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and 

A12 = 2uC, (26d) 

B l2 = -2vC, (26e) 

C = -u(A 12 -B 12 ), (26f) 

with the only nonvanishing initial conditions ^4(0) = fz^p + n th and ^12(0) = j~^t- The 
solutions now read 

,,2 



= [A 2 + n,,(l-A 2 )]^, (27a) 

B(t) = [X 2 + n th (l-X 2 )]f^, (27b) 

D(t) = \[X 2 + n th {l-X 2 )]^^-, (27c) 

, . . , cos 2 vt 

A 12 (t) = X— 

„ . , sin 2 vt 

B l2 (t) = X^—^, 

_ -A sin(2z/t) 

C[t) ~ ~~2 T^X?~' 



(27d) 
(27e) 
(27f) 
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C. Fidelities - Parametric Interaction 



Using the results of Appendix the integrals (fTij) can be easily evaluated giving the 
following results 



7(0:1,02; 0, 0) 
/(oi,o 2 ;0, 1) 



exp [—A (|oi| 2 + |o 2 | 2 ) + Au (oio 2 + 0^02)] 



(28a) 



B + l 



exp 



A 



D 2 
B + l 



Ol 



C 1 



M2 + -77-7-7 ) (oio 2 + ajog) 



B + l 



I(oi,a 2 ; 1,0) 



1 



B + l 



exp 



.4 



C 1 



B + l 



Ol 



A 



B + l 



D 2 
B + l 



Q 2 



(28b) 



Q 2 



+ [A 12 + 



CD 
B + l 



(0102 + 0JO2) 



(28c) 



1(01,02; 1,1) 



(B + I) 2 - B 2 



■ exp 



12 



[B + l){C 2 + D 2 ) + 2B 12 CD \ , 2 2 ^ 
A — — ^ I (|oi| + |o 2 | J 



(B + l) 2 -B 2 



12 



2{B + l)CD + B 12 (C 2 + D 2 ) . 
+ I A12 H 77; ; 773 ) (0102 + aia 2 ) 



(B + I) 2 - B\ 



12 



(28d) 



The probabilities are then derived using the condition ([13 



Poo 
P01 
P10 
Pll 



1(0,0; 1,1) 



(B + I) 2 - B\ 2 ' 



1(0,0; 1,0) -1(0,0; 1,1) 
1(0,0; 0,1) -1(0,0; 1,1) 



1 



B + l (B + l) 2 -B 2 2 
1 1 



B + l (B + l) 2 -B 2 2 
1(0, 0; 0, 0) - 1(0, 0; 1, 0) - 1(0, 0; 0, 1) + 7(0, 0; 1, 1) 
1 



1 



+ 



B + l (B + l) 2 -B 2 



(29a) 
(29b) 
(29c) 



(29d) 



12 



Finally, with the help of Eqs. (ll2p and (I18p we can easily arrive at Eq. (jl9j) . 
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D. Fidelities - Beam Splitter Interaction 



Using the results of Appendix [Bj the integrals (|14p can be easily evaluated giving the 
following result 



7(0:1,02; 0, 0) 
7(oi,o 2 ;0, 1) 



exp [—A (|oi| 2 + |o 2 | 2 ) + Au (oio 2 + 0^02)] 



(30a) 



B + l 



exp 



A 



C 2 



B + l 



Ol 



D 2 
B + l 



Q 2 



^12 + 77——, ) (oio 2 + 0^03) 



B + l 



(30b) 



1(01,02; 1,0) 



1 



B + l 



exp 



A 



D 2 
B + l 



■i\ 2 ~ i A --^—r ) l«2| 2 



B + l 



+ [A-12 + 



CD 
B + l 



(0102 + 0JO2) 



1(01,02; 1,1) 



■ exp 



.4 



{B + l){D 2 +C 2 ) + 2B 12 CD\ 



(B + I) 2 - B 2 2 

2(B + l)CD + B 12 (C 2 + D 2 



(B + I) 2 - B\ 2 



J (K| 2 + |a 2 | 2 ) 



+ [A 12 + 



(B + I) 2 - B\ 



12 



(0102 + 0^02) 



(30c) 



The probabilities are then derived using the condition ([13 



poo = 7(0,0; 1,1) 



P01 
P10 



(B + I) 2 - B 2 2 ' 



7(0,0; 1,0) -7(0,0; 1,1) 
7(0,0; 0,1) -7(0,0; 1,1) 



1 



B + l (B + l) 2 -B 2 2 
1 1 



B + l (B + l) 2 -B 2 2 
7(0, 0; 0, 0) - 7(0, 0; 1, 0) - 7(0, 0; 0, 1) + 7(0, 0; 1, 1) 
1 



1 



+ 



B + l (B + l) 2 -B 2 



(31a) 
(31b) 
(31c) 



(31d) 



12 



Finally, with the help of Eqs. (ll2p and (I18p we can easily arrive at Eq. (jl9j) . 
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Figure Captions 

Figure 1: Two-entangled light fields, generated for instance by an Optical Parametric 
Oscillator (OPO), impinge each on one atomic ensemble. The light field couples to one 
dipolar transition of a A-configuration of electronic levels (coupling strength g), while the 
second transition is driven by a laser (coupling strength f2). The light fields are here labeled 
by the corresponding photon-annihilation operators a\ and 02, while the involved atomic 
levels are \g) and |e) (stable states), with \i) the common excited state. Before the interaction 
the atoms are prepared in state \g). The effective interaction between fields and collective 
atomic excitations is parametric for the coupling scheme shown in (a), namely, photonic 
and atomic excitations are simultaneously created or annihilated. In (b) the interaction is 
beam-splitter-like, i.e., the total number of photonic and atomic excitations is conserved. 
After the interaction, the collective state of the atomic ensembles is measured. As a result 
the entanglement between fields a\ and 02 is expected to increase (see text for details). 

Figure 2: Fidelities for the four conditional states, obtained after parametric interaction 
with the atomic ensemble from Eqs. f)19|> . versus the interaction time t (in units of l//x). 
Here, A = 0.5 and nth — 0. The curves correspond to Foo (dot-dashed line), Fn = F\q 
(dashed line) and F\\ (solid line). The initial teleportation fidelity i^ n it (dotted line) is 
shown for comparison. 

Figure 3: Fidelities for the beam-splitter interaction, as in Eqs. f)19[) . as a function of 
A, and for various values of nth- The fidelities are evaluated at the optimal time t ~ vr/z^, 
where Fn is maximal. The curves correspond to Foo (dot-dashed line), Fqi = F\q (dashed 
line) and Fu (solid line). The initial teleportation fidelity Fj n jj (dotted line) is shown for 
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comparison. 

Figure 4: Efficiency £, Eq. (I21j) as a function of t (in units of and A for the beam- 
splitter interaction and (a) n t h = 0, (6) = 0.05. 

Figure 5: Fidelities for entanglement swapping using the parametric interaction as a 
function of fit for various values of A and n^. The curves correspond to Fqq (dot-dashed 
line), Fqi = F\q (dashed line) and F\\ (solid line). The initial teleportation fidelity i^nit 
(dotted line) is shown for comparison. 

Figure 6: Fidelity for entanglement swapping using the parametric interaction as a 
function of fit for various values of A. This fidelity is compared with the entanglement 
stemming from a pure parametric interaction, F param = (1 + fit)/2. 
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Fig. 1. Two-entangled light fields, generated for instance by an Optical Parametric 
Oscillator (OPO), impinge each on one atomic ensemble. The light field couples to 
one dipolar transition of a A-configuration of electronic levels (coupling strength 
g), while the second transition is driven by a laser (coupling strength Q). The light 
fields are here labeled by the corresponding photon-annihilation operators a± and 
02, while the involved atomic levels are \g) and |e) (stable states), with \i) the 
common excited state. Before the interaction the atoms are prepared in state \g). The 
effective interaction between fields and collective atomic excitations is parametric 
for the coupling scheme shown in (a), namely, photonic and atomic excitations are 
simultaneously created or annihilated. In (b) the interaction is beam-splitter-like, 
i.e., the total number of photonic and atomic excitations is conserved. After the 
interaction, the collective state of the atomic ensembles is measured. As a result the 
entanglement between fields a\ and 02 is expected to increase (see text for details). 
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Fig. 2. Fidelities for the four conditional states, obtained after parametric interaction 
with the atomic ensemble from Eqs. (|19p . versus the interaction time t (in units of 
l//i). Here, A = 0.5 and n t h = 0. The curves correspond to Fqq (dot-dashed line), 
Fq\ = Fio (dashed line) and Fu (solid line). The initial teleportation fidelity Fi n a 
(dotted line) is shown for comparison. 
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Fig. 3. Fidelities for the beam-splitter interaction, as in Eqs. (fl~9l) . as a function of 
A, and for various values of nth- The fidelities are evaluated at the optimal time 
t ~ vr/z^, where F\\ is maximal. The curves correspond to Fqq (dot-dashed line), 
Foi = Fio (dashed line) and F\\ (solid line). The initial teleportation fidelity Fi n a 
(dotted line) is shown for comparison. 
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(b) 




4. Efficiency S , Eq. (|2ip as a function of £ (in units 
l-splitter interaction and (a) = 0, (b) n t h = 0.05. 
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Fig. 5. Fidelities for entanglement swapping using the parametric interaction as a 
function of \xt for various values of A and nth- The curves correspond to Foo (dot- 
dashed line), Fqi = F±q (dashed line) and Fix (solid line). The initial teleportation 
fidelity i^ n & (dotted line) is shown for comparison. 
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X= 0.01 




Fig. 6. Fidelity for entanglement swapping using the parametric interaction as a 
function of fit for various values of A. This fidelity is compared with the entanglement 
stemming from a pure parametric interaction, F param = (1 + /it)/2. 
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